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ABSTRACT 


A  navigation  scheme  is  investigated  that  uses  a  Schuler- tuned  inertial  platform  that  is 
maintained  orthogonal  to  the  mass-attraction  gravity  vector.  The  expressions  for  platform 
acceleration  and  velocity  are  developed,  and  a  mechanization  to  achieve  navigation  is 
determined.  Results  indicate  that  an  inertial  platform  that  is  not  torqued  in  the  vertical 
axis  can  be  readily  mechanized  to  seek  a  mass-attraction  vertical  and  to  accomplish 
navigation  on  earth. 
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Publication  of  this  technical  documentary  report  does  not  constitute  Air  Force  approval 
of  the  report’s  findings  or  conclusions.  It  is  published  only  for  the  exchange  and  stimula¬ 
tion  of  ideas. 
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INTRODUCTION 

To  navigate  in  the  geographic  coordinate  system'  of  the  earth,  we  desire  to  use  acceler¬ 
ometers  that  are  mounted  on  an  inertial  platform.  The  total  acceleration  in  inertial  space 
as  experienced  by  the  platform  that  is  positioned  at  a  point  P.  which  is  a  distance  R  from 
the  center  of  the  earth,  will  be  derived.  A  particular  orientation  of  the  platform  is  chosen 
so  that  its  vertical  axis  is  coincident  with  the  vector  ^  from  the  center  of  the  earth.  In 
the  development,  the  geographic  center  of  the  earth  and  the  center  of  mass  attraction  will 
coincide  so  that  the  platform  will  be  aligned  to  the  mass-attraction  gravity  vector. 

The  position  of  the  platform  on  the  earth  will  be  expressed  in  terms  of  the  geographic 
longitude  angle  and  the  mass-attraction  latitude  angle.  Note  that  in  contrast  the  geographic 
latitude  is  measured  from  a  plumb  bob  vertical.  The  mass  suspended  by  a  plumb  line  is 
attracted  by  the  gravitation  of  the  earth  and  is  influenced  by  the  rotation  of  the  earth  to 
cause  the  plumb  line  to  point  along  the  vector  sum  of  the  mass-attraction  vector  and  the 
centripetal  acceleration  vector.  To  account  for  a  discrepancy  in  latitude  position  resulting 
from  a  determination  of  position  using  mass- attraction  vertical  and  a  determination  of 
position  using  a  plumb  bob  vertical,  we  will  alter  the  geographic  inputs  and  outputs  to  the 
platform  for  the  contribution  of  centripetal  acceleration. 

To  mechanize  the  expressions  derived  for  acceleration  and  velocity  of  the  platform,  we 
will  use  the  Schuler  principle  of  developing  an  angular  rate  of  the  platform  that  is  equal 
to  the  inertial  velocity  divided  by  the  radius  of  the  earth  to  maintain  a  mass-attraction 
vertical  and  to  navigate  over  the  earth. 


DERIVATICa'J  OF  MASS  ACCELERATION 
IN  A  ROTATING  COORDINATE  SYSTEM 

_  A  right-hand  coordinate  system  is  selected  with  Us  or^in  at  the  center  of  the  earth,  its 
K  axis  through  the  North  Pole  of  the  earth,  and  its  I  and  J  axes  in  the  plane  of  the  earth’s 
equator  with  I  passing  through  the  Greenwich  meridian.  In  this  development,  the  radius  of 
the  earth  is  considered  constant  for  all  angles  above  or  below  the  equatorial  plane.  See 
Figure  1. 


1  NORTH 

K 

X 

/ 

J 

/■ 

I 

EQUATOR 

Figure  1.  Coordinate  System  With  Origin  at  Center  of  Earth 
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We  desire  to  determine  the  acceleratmn  in  inertial  space  that  is  experienced  at  point  P 
while  P  is  free  to  move  in  the  I ,  J,_and  K  coordinate  system  and  while  the  1,  J,  and  K 
axes  rotate  in  inertial  space  about  K  with  an  angular  rate  cog.  In  the  development,  the 
following  terms  will  be  used: 


w  K 
e  e 

R  - iR  +  TR  +  KR^ 

X  y  z 

R  ■  inertial  velocity  of  P 

Vp  ■  velocity  of  P  relative  to  1,  J.  and  K  axes 

R  ■  inertial  acceleration  of  P 

ap  ■  acceleration  of  P  relative  to  1,  J,  and  K  axes 

Aj  ■  acceleration  of  P  in  inertial  space  including  gravity  of  mass  attraction 

g  ■  gravity  vector  of  mass  attraction  of  earth  assumed  directed  toward  the 
center  of  the  earth. 

will  be  differentiated  once  to  form  R,  the  inertial  velocity  of  P,  and  once  again  to  form 
,  the  inertial  acceleration  of  P.  Therefore. 


R 

-  IR  +  JR 

+  KR 

X  y 

z 

R 

-  iR  +  JR 

4  KR 

X  y 

z 

R 

-  w  X  R  + 

V 

e 

P 

R 

■  U  X  R  + 

P 

Wp  X  i 

KR 


R 

R 


V  -  IR  +  JR  +  KR  +  IR  +  JR  +  KR 
p  X  y  z  X  y  z 

_ 

V  «  ui  X  V  +  a 
p  e  p  p 

^  «  O.  CO  considered  constant 
e  e 


CO  X  R  +  V 

e  p 

CO  x(co  xR+V)+co  xV  +a 

e  '  e  P'  e  P  P 


R  ■0)  xco  xR  +  2cl'  X  V  +  a 
e  e  ®  P  P 


A,  -  X  co„  X  R  +  2iji  X  V„  +  a„  +  g„ 
lee  e  p  p  '’m 
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The  preceding  expression  for  A.  is  the  acceleration  of  P  in  inertial  space,  which  differs 
from  R  by  the  gravity  term  g^.  A  plumb  bob  on  the  earth  will  be  under  the  influence  of 
A,  to  establish  a  vertical  direction,  which  is  the  vector  sum  of  R  +  g  .  Considerations 
of  the  acceleration  of  center  of  the  earth  relative  to  the  center  of  inertial  space  are 
neglected  in  A^  as  insignificant. 

A-  -  g  in  terms  of  I,  J,  and  K  is  as  follows: 

I  m 

CO  xR«I(-coR)+J(coR) 

0  \  0  y/  '  e 

CO  xw  xR=1(-co®R)+J(-co®R) 
e  e  '  e  'ey' 

2w  X  V  »  2l  (-to  R  )  +  2J  (to  R  ) 
e  p  '  e  y'  '  e  x^ 

A  -  1  r  (-to-  ®  R  )  -^  2(-to  R  )  R  1 
1  L '  e  x'  'ey'  x  J 

*  j  [(-W  *  *  "y  ] 

*K  [RJ 

m 

ORIENTATION  OF  A  COORDINATE  SYSTEM 
ALONG  THE  MASS-ATTRACTION  VECTOR 

We  desire  to  orient  a  coordinate  system  i,  j,  and  k  at  P  so  that  k  is  along  R,  i  is 
directed  toward  K  (North),  and  j  then  is  directed  eastward.  See  Figure  2. 


Figure  2«  Coordinate  System  With  One  Axis  Directed  Toward  Center  of  Earth 
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In  this  development,  the  following  terms  will  be  used: 

L  ■  angle  eastward  from  Greenwich,  longitude 
X  ■  angle  from  the  equatorial  plane,  mass-attraction  latitude 
ic  ■  1  cos  X  cosL  +  J  cos  X  sinL  +  K  sin  X 
J  ■!  (-SinL)  +  J  (cosL) 

i  ■!  (-sin  X  cosL)  +  J  (-sin  X  sinL)  +  K  cos  X. 

The  intent  of  positioning  k  to  be  coincident  with  R  is  to  provide  a  coordinate  system 
with  one  axis  directed  toward  the  center  of  the  earth.  When  the  geographical  center  of 
the  earth  is  assumed  to  be  located  at  the  center  of  mass  attraction,  i,  j,  and  k  will  then 
be  aligned  along  the  mass-attraction  vector. 

The  acceleration  of  point  P  is  now  to  be  expressed  in  terms  of  along  1,  j,  and  k 
as  follows: 


I  [(<»>*  R  sin  X  cos  L  +  a-  ’  R  sin  X  sin  L) 

L'  e  X  ^  y  ^ 

+  (2a-  R  sin X  cosL  -  2a:  R  sinX  sinL) 
'ey  ex  ' 

••  ••  •• 

(R  sinX  cosL  R  sinX  sinL  -  R  cosX) 
'  X  y  z  ' 


j  f  (a*  *  R  sin  L  -  cc  ’  R  cos  L) 
L'  e  X  e  y  ' 

+  (2a:  R  sinL  +  2a-  R  cosL) 
'ey  ex  ' 

••  •• 

+  (R^  sin  L  -  Ry  cos  L)  j 


+  k  [-(a;  ®R  cosXcosL  +  «  ®R  cos X  sinL) 
L  '  e  X  e  y  ' 

•  • 

-  (2u:  R  cos  X  cos  L  -  2a.'  R  cos  X  sin  L) 
'ey  ex 


To  express  Aj  in  terms  of  R,  R,  and  R,  we  must  derive  the  following  equalities: 


4 


ASD-TDR-63-168 


*  RcosA  cosL 
Ry  ■  RcosX  sinL 
R^  «  RsinX 

R  cosL  *  R  SinL 

y  .  X 

R  cosL  -  R  sinL  *  0 

y  X 

R  cosL  +  R  SinL  «  RcosX 

X  y 

R  cosL  +  R  SinL  =  RcosX  -  RX  sinX 
X  y 

R  cosL  -  R  SinL  *  RLcosX 

y  X 

(R^cosL  RySinL)sinX  *  -  RsinX  cosX  [(L)®  ^ 

-  sin®  X  (2RX  +  RX)  R  cosX  sinX 

••  ••  m  9 

(R^cosL  +  RySinL)cosX  =  -  Rcos®  X  [(L)®  +  (X)®] 

-  cobX  SinX  (2RX  +  R  X )  -^  R  cos®  X 

•  •  •«  ••  •  •  •  ■ 

R  cosL  -  R  SinL  =  RL  cosX  -  2RLX  sinX  2RLcosX 
y  ^ 

R^  =  RsinX  +  RX  cosX 

R^  =  R  SinX  *  2RX  cosX  +  RX  cosX  -  R(X)®  sinX 
Then  substituting  the  equalities  gives 


1  [RsinX  cosX  (0;^  +  L)®  +  2RX  +  RX  ] 

+  j  [  2(u.‘g  +  L)  RcosX  -  2(u!g  +  L)RX8inX  +  RLcosx] 
+  k  [-(a!g  +  L)®Rco8®X  -R(X)*  +  «  +  ] 
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MECHANIZATIC*^  OF  ACCEL  ERA  ri(»J  EXPRESSIONS 
TO  ACHIEVE  INERTIAL  VELOCITY  OF  POINT  P 

We  desire  to  integrate  the  acceleration  expressions  derived  so  that  the  velocity  of  P  in 
inertial  space  is  expressed  in  a  coordinate  system  Ti  ,r^  ,  and  that  is  related  to  i,7f 
and  k  as  follows: 

r,  «  k 

■  i  sin  ^  +  j  cos  ^ 

P  P 

■  i  cos  ^  -  j  sin  ^  , 

P  P 


ra 


Figure  3  shows  the  desired  coordinate  system. 


Figure  3.  Coordinate  System  r.,  r.,,  and  r 

Angle  measured  counterclockwise  from  i  (North),  at  this  step  in  the  development  will 
not  be  restricted  to  a  magnitude  or  rate  ('^p).  Eventually  ,  and  will  be  restricted 

so  that  there  will  be  no  rotation  of  the  coordinate  system  about  r^  in  inertial  space. 


Aj  is  now  expressed  in  the  r^ ,  r^  ,  and  ra  directions: 


ri  [-(a'g  +  L)®Rc08®X  -R(X)®  +  R 


+  r,  f  (u  +  L)’ Rsin  X  cos  X  sin 'I'  +  2RX  sin'l' 

^  L  '  e  '  p  p 

+  2(0.’  +  L)  Rcos  X  cos  ^  +  RX  sin  ^ 

'  e  '  p  p 

-  2(o.’g  +  L)RX  sinX  cos'I'p  +  RLcosX  cosi'p  ] 
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—  r  *  *  ’ 

*  (<«■•  +  L)*  Rsin  X  cos  X  cos  'I'  +  2RX  cos  ^ 

L  '  e  '  p  p 

•  •  •  • 

-  2(cu  +  L)  Rcos  X  sin  ^  +  RX  cos 
e  P  P 

•  •  ••  ^ 

2(a,'  +  L) R X  sinX  sin'*'  -RLcosXsin^'  I 
'  e  '  p  p  •* 

Upon  integration  of  Aj,  the  inertial  velocity  of  P  should  be  equal  to  which  was  pre¬ 
viously  derived.  Therefore, 

R  ■  w  X  R  V 
e  p 

X  R  ■  1  (-‘^'gRy)  +  J 

u)  X  "r  «  j  u'  RcosX 

V  -iR  +  JR  +KR 

p  X  y  z 

V  ■  i  RX  +  j  RLcosX  +  k  R 
P 

R  ■  i  R  X  +1  (Wg  +  L)  Rcos  X  +  k  R 

R  ■  Ta  [  RXcos^'p  -(a,'g  +  L)  RcosX  sin’I'p  j 
H-Tj  [RXsin'i'  +  (a-  +  L)  Rcos  X  cos 'i'  ] 

+  [r] 

R  is  to  be  the  output  of  three  isolated  integrators  following  three  accelerometers  aligned 
along  r, ,  Tg  ,  and  ,  respectively,  which  sense  Ap  Differentiating  R  to  achieve  R  will 
give  the  expression  for  the  accelerations  required  at  the  input  to  the  integrators.  Therefore, 

1  - 

ta  [rXcos'I'p  -(w^  L) RcosX  Sin ifp  ] 

[•  •  ••  •  • 

RX  cos  '*'  +  RX  cos  ^  -(u3  +  L)  Rcos  X  sin  ♦ 
p  p  '  e  '  p 

•  •  •  • 

*  (w  +  L)  RX  sin  X  sin  'I'  -  L  Rcos  X  sin  'I' 

'  e  '  p  p 

•  •  •  •  ^ 

-  RX  -((c»  +  L)R^  cosXcos^  I 

P  p  '  e  '  p  p 
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+  r j  ^  RX  sin  'i'p  +  +  L)  Rcos  X  cos  '1'^  ] 

—  r  *  *  *  *  • 

+  tj  RX  sin^  +  RX  sin'J'  +  (a,-  +  L)RcosX  cos^' 

L  p  p  '  e  '  p 

•  •  •  • 

“  (o:  +  L)  RX  sin  X  cos  ^  +  L  Rcos  X  cos  ^ 

^  P  P 

•  •  •  • 

RX 'I'  cosv|>  -(u;  +  L)R'*'  cos  X  sin  1 

P  p  '  e  ^  p  p  J 

*j.  [r]  *?,  [ii] 


in  the  expression  for  R  is  evaluated  as  follows: 


Cl'  *  co^K 
e  e 


a'  *  i  a*  cosX  +  k  a’  sinX 
e  e  e 


From  the  preceding  expression,  a  rotation  a*^sinX  occurs  about  k  in  inertial  space.  The 


restriction  that  the  rj  ,  Tg  ,  and  r,  coordinate  system  cannot  rotate  in  inertial  space  about 
r^  will  now  be  incorporated  i 
the  angle  between  i  and  r,  , 


r^  will  now  be  incorporated  in  the  development.  When  r^  is  coincident  with  k  and  is 


^  =  ~  (a'  +  L)  sinX  . 

p  '  e  ' 

The  rate  of  is  negative  because  positive  rotation  was  selected  counterclockwise  for 

To  accomplish  no  rotation  in  inertial  space  about  r,  ,  we  will  not  torque  the  inertial  plat¬ 
form  in  this  axis. 

The  r  terms  in  R  are  evaluated  as  follows: 

oJ  R  =  r«,  RXcos^  -  (cc  L)  R  cos  X  sin  >1^  J 
P  ^  ^  P  ^  ^  '  P^ 


+  r 


r^  [  RXsin^  ^  L)  R  cos  X  cos ] 

F.  [n 


where  cl'^  is  the  rotation  of  r^,  r^,  and  r3  in  inertial  space.  To  accomplish  navigation  in 

the  coordinate  system  of  the  earth,  we  will  follow  the  Schuler  principle  for  generation  of 

Wp.  In  a  Schuler  mechanization,  the  integrated  acceleration,  which  is  velocity,  for  each  of 

the  two  horizontal  axes  (r,  and  r-,)  is  divided  by  R,  the  radius  of  the  earth,  which  gives 
an  angular  rate,  cc,  for  torquing  the  inertial  platform  to  maintain  the  aforementioned  axes 
in  the  horizontal  plane.  Following  this  procedure,  we  obtain 
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"2  '2  '  -  -r  '3 


Note  that  a.*,  is  negative  to  maintain  r,  horizontal. 

The  Vj,  Vp,  and  V,  inertial  velocities  are  found  from 

V,  ■  R  ^  X  sin  +  (Wg  +  L)  cos  X  cos  '1'^  j 

V,  ■  R  [  X  cos 'I'  -(a;  +  L)  cosX  sin 'I'  ] 

V,  .  [b] 

Note  that  ojiTj  is  zero,  therefore 
OJp  •  03s  Ta  UJ,  r, 

i  -  WsTs  *  v,7,  *  W,r, 

a’p  X  ^  ^  a;*  Vj  -  <i.'3Ve|  rj  +  ' 

cUp  X  R  ■  Ti  [-R(X)®  -(ct'^  +  L)*Rcos*  x] 

*  [(«)]  [  X  sin'l'p  +  (a'g  +  L)  cos  X  cos 'J'p  ] 

*  r,  [(R)  ]  [x  cos'l'p  -(a'g  +  L)cosX  sin'J'p  ] 

Substituting  expressions  developed  for  ^p  and  r  terms  gives 

••  r  •  *  **  *  * 

R  ■  rj  [  RX  cos 'I'  +  RX  cos 'I'  -(w  +  L)RcosX  sin 'I' 

P  p  6  P 

•  •  •  • 

+  2(u'  L)  R  X  sin  X  sin  ^  -  L  Rcos  X  sin  '!> 

'  e  '  p  p 

+  (Wg  +  L)*  Rsin  X  cosX  cos'l'p  ] 

+  r,  f  RX  sin 'I'  *  RX  sin’p  +  (cc  *  L)RcosX  cos^ 

»i  p  p  ''  e  '  p 

•  •  •  • 

-2(o3  +  L)RX  sinX  cos'l'  +  L  RcosX  cos’l' 
e  '  p  p 

+  (cc  +  L)®  RsinX  COSX  sin 'I' „  J 
'  e  '  p 

-.[«i 


+  cc„  X  R 
P 
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If  the  value  of  Wp  x  R  were  included  in  the  Tj,  r^,  and  r,  directions  above  for  R,  we 


would  see  that  R 


1  m 


At  this  point  in  the  development,  we  must  examine  Aj  and  the  accelerometers  in  the 
rj,  Tg,  and  r,  coordinate  system.  For  accomplishment  of  the  Schuler  mechanization,  the 
platform  is  rotated  in  inertial  space  by  Wp.  Torque  is  applied  to  the  platform  at  a  rate 
u'p  in  the  two  horizontal  axes  to  oppose  the  inertial  position  sought  by  the  platform  so 
that  an  earth  orientation  can  be  maintained.  Because  the  platform  and  the  accelerometers 
have  an  inertial  velocity  R,  a  rotation  of  the  platform  coordinate  system  in  inertial  space 
results  in  an  apparent  acceleration  of  the  platform  equal  to  a^p  x  R  as  seen  in  the  coor¬ 
dinate  system  of  the  earth.  The  contribution  of  uj  x  R  to  R  is  determined  from  the  total 
_  _  P 

integral  of  Aj  -  g^: 


R  *  J*(  Aj  - 


■[/ 


(A,  -  gm>  I*  n  + 


]  7,  +  [/A2<lt]r2 


[/A3df]73  + 


J twpX  R)  df 


Note  that  the  above  expression  for  Aj  was  verified  by  differentiating  R, 

The  preceding  expression  clearly  shows  the  components  of  the  total  Aj  that  will  be 
sensed  by  accelerometers  aligned  along  r, ,  r^,  and  r^,  namely  A^,  Ap,  and  A*,  respec¬ 
tively.  To  accomplish  the  total  integration  of  A^,  we  must  supply  the  x  ^  term  at 
the  input  to  the  three  integrators  receiving  A^,  Ap,  and  In  the  mechanization  of  the 

desired  inertial  navigator,  therefore,  the  velocity  of  point  P  can  be  achieved  from  the 
accelerometers  by  adding  the  following  quantities: 

To  accelerometer  output  Ap  must  be  added 

[(R)]  [xcosO'p  -(u'g  L)cosXsinM'p  ] 

which  is  equal  to 

To  accelerometer  output  Ag  must  be  added 

[TO]  [  X  sin4j^  +  (cc^  +  L)cosX  cos^I^p  ] 

which  is  equal  to 
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To  accelerometer  output  A,  must  be  added 

[-R(X)>  -(o-'g  +  L)®Rcos®  X  ] 
which  is  equal  to  w-Va  -  o-'aV-  -  g„. 

In  the  closed-loop  operation  of  the  inertial  platform.  Figure  4,  the  a’p  and  R  terms  that 
must  supply  a;  x  R  are  readily  available  as  outputs  of  the  platform. 


DETERMINATION  OF  PLATFORM  POSITION  IN  THE 
COORDINATE  SYSTEM  OF  THE  EARTH 


Navigation  in  the  coordinate  system  of  the  earth  is  accomplished  by  using  the  torquing 
rates,  Wj  and  u-'a,  that  are  generated  in  the  mechanization  of  the  Schuler- tuned  inertial 
navigation  system  as  follows; 

Mass  attraction  latitude  results  from  an  integration  of  X ,  where. 


X  ■  tL'a  cos  +  cc,  ein  , 
P  P 

and  longitude  results  from  an  integration  of  L,  where. 


ta.',  cos  'I'  -  CO-  sin  'I'  -  w 
_ p  P  ^ 

cosX 


In  these  expressions,  co  ,  rotation  of  the  earth,  is  a  constant,  and  'T”  is  determined  from 
.  e  p 

an  integration  of  ^p,  where. 


P 


-(cOg  +  L)sinX. 


Mass  attraction  latitude  is  transformed  to  geographic  latitude  as  follows: 


g 


a  + 
c 


’m 


where 


g  ■  plumb  bob  vertical 

a^  -  earth’s  centripetal  acceleration  component  that  is  perpendicular  to  g 

i  -  mass- attraction  vertical 
®m 

a  ■  angle  between  g  and  g^^. 


m 


Note  that  a<  1  ,  then 


( ***>  )*  R  »in  X  cot  X 


(N  RADIANS. 


“  flm 

Angle  a  is  added  to  the  mass-attraction  vertical  angle  for  latitude  to  obtain  geographic 
latitude. 
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Figure  4,  Mechanization  of  Velocity  Expressions 
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CONCLUSIONS 

The  mechanization  of  inertially  measured  accelerations  to  derive  inertial  velocity  and 
to  accomplish  navigation  on  earth  is  readily  achieved  from  three  accelerometers  that  are 
orthogonally  mounted  in  a  platform  coordinated  system,  which  is  vertically  aligned  to  the 
mass-attraction  vector  of  the  earth  and  which  is  not  torqued  about  the  vertical  axis. 

To  obtain  platform  velocity  relative  to  the  earth,  we  must  subtract  the  velocity  of  the 

earth,  x  R,  from  the  inertial  velocity  output  of  the  platform.  Because  latitude  position 

that  is  determined  by  the  inertial  navigator  does  not  include  centripetal  acceleration,  an 
angular  displacement  must  be  added  to  the  latitude  output  of  the  platform  to  obtain  geo¬ 
graphic  latitude. 

The  radius  of  the  earth  was  assumed  constant  in  the  development.  Varying  R  as  a  func¬ 
tion  of  latitude  can  account  for  the  spheroidal  shape  of  the  earth.  To  include  the  altitude 
of  the  platform  above  the  surface  of  the  earth  in  the  total  distance  of  the  platform  from 
the  center  of  the  earth,  we  can  add  an  integration  of  platform  velocity  in  the  vertical, 

ri,  direction  to  R. 
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